Observation of a multiferroic critical end point 
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The study of abrupt increases in magnetization 
with magnetic held known as metamagnetic tran- 
sitions has opened a rich vein of new physics in 
itinerant electron systems, including the discov- 
ery of quantum critical end points with a marked 
propensity to develop new kinds of order. How- 
ever, the electric analogue of the metamagnetic 
critical end point, a "metaelectric" critical end 
point has not yet been realized. Multiferroic ma- 
terials wherein magnetism and ferroelectricity are 
cross-coupled are ideal candidates for the explo- 
ration of this novel possibility using magnetic- 
field (H) as a tuning parameter. Herein, we 
report the discovery of a magnetic-field-induced 
metaelectric transition in multiferroic BiM^Os 
in which the electric polarization (P) switches 
polarity along with a concomitant Mn spin-flop 
transition at a critical magnetic field H c . The si- 
multaneous metaelectric and spin-flop transitions 
become sharper upon cooling, but remain a con- 
tinuous crossover even down to 0.5 K. Near the 
P=0 line realized at // #c«18 T below 20 K, the 
dielectric constant (e) increases significantly over 
wide field- and temperature ( T)-ranges. Fur- 
thermore, a characteristic power-law behavior is 
found in the P(H) and e(H) curves at T=0.66 
K. These findings indicate that a magnetic-field- 
induced metaelectric critical end point is realized 
in BiM^Os near zero temperature. 

The term "critical end point" refers to a singular point 
in the phase diagram of matter at the end of a 1st order 
phase line, as for example, the liquid-gas critical point 
of water. The importance of this special point for broad 
classes of matter has rapidly increased in recent years 
[l], Not only can it provide large thermal fluctua- 
tions as a necessary ingredient for displaying universal 
power-law of physical quantities, but in the special case 
where the phase transition is suppressed to zero tem- 
perature, it can give rise to intense quantum mechanical 
fluctuations with a marked propensity to develop insta- 
bilities into new kinds of ground state. The latter case 
has been recently realized in itinerant metamagnets such 
as Sr 3 Ru2 07 [3| and URu2Si2 [J]. These developments 
motivate a parallel search for a metaelectric critical end 
point. By analogy with magnetism, one might expect 



such sudden changes in electric polarization to be a rather 
general phenomenon in (anti)ferroelectric materials un- 
der application of an electric field (E). However, to date, 
those effects have been limitedly observed in specific sys- 
tems such as relaxor ferroelectrics [5[ and DyVC>4 with 
the Jahn- Teller structural distortion [6]. One possible 
reason for the scarcity of the phenomenon is the practi- 
cal difficulty of applying the large voltages required (^1 
kV) without inducing electrical breakdown. Magnetic 
fields may in fact be better candidates for inducing meta- 
electric transitions, for not only do they avoid the prob- 
lem of electrical breakdown, they provide a reversible 
method of tuning matter with many well-established ad- 
vantages over alternative methods such as pressure or 
chemical doping. This is still a challenging task be- 
cause of the intrinsically small cross-coupling between 
spin and lattice degrees of freedom. However, in a spe- 
cial class of matter called multiferroics where the cross- 
coupling between electricity and magnetism is enhanced, 
0, S, H, EL E3 ! such a rare metaelectric transition, 
and as its critical end point can be induced by application 
of magnetic field, as we show in this study. 

Among recently discovered multiferroic materials, the 
family of i?Mn20s (i?=rare earth, Y) are of particular 
interest, because of their large dielectric response under 
magnetic fields Q. While B1MTI2O5 is isostructural to 
i?Mn205 (Fig. la), its magnetic ground state is dis- 
tinct from those of other fflVh^Os compounds. Only 
BiMn205 remains in a commensurate antiferromagnetic 
(AF) phase below Tn~40K with a propagation vector of 
Q=(0.5, 0, 0.5) [HGj], while other -RMn 2 5 compounds 
develop an instability into an incommensurate (IC) AF 
structure below Xic«25 K[l3j]. According to DC magne- 
tization (M) (Fig. 2a) and neutron diffraction data [f4| . 
BiMn205 has a typical AF spin ordering with an easy 
axis almost parallel to the a-axis without showing any 
signature of a spiral ordering. At Tn 7 the development 
of magnetic order is accompanied by the growth of fer- 
roelectric (FE) polarization P and a sharp peak in the 
dielectric constant e (Figs. 2b and 2c). 

Multiferroicity in materials with incommensurate spi- 
ral spin ordering is normally accounted for in terms of 
the "spin-current model" [151 ] . Recent neutron scatte ring 
studies of i?Mn 2 5 (R=Tb, Ho, Y, and Bi) 0, []J, [if 
have suggested an alternative "exchange-striction" see- 
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FIG. 1: Magnetic structure and phase diagram, (a) The 
structure and AF zigzag spin chains under various magnetic 
fields for BiMn2C>5. J 3 and J 4 denote the interchain and 
intrachain exchange interactions between Mn 3+ -Mn 4+ spin 
pairs, respectively, while J 5 describes that between neighbor- 
ing Mn 3+ spins in the bi-pyramid. Red, dark blue, and light 
blue arrows represent spin directions at H—0, H=H C , and 
H^H C , respectively. The corresponding simplified spin ar- 
rangements and the net P (black open arrows) are shown for 
(6) H=0, (c) H=H C , and (d) H^>H C . (e) The H-T phase 
diagram and the intensity contour of P based on the P{H) 
data. All symbols were determined from the peak positions of 
e(H), e(H), and dM/dH. (/) Suggested schematic phase dia- 
gram. Dotted, solid, and dashed lines indicate the 2nd-order 
paraelectric (PE)-FE transition, lst-order AF to spin-flopped 
AF (AF-SF) transition, and AF to AF-SF crossover, respec- 
tively. A possible critical end point T* is denoted as a solid 
circle. 



nario for this commensurate magnetic system. The AF 
phase of -RM^Os involves zigzag chains of spins (dashed 
orange lines in Fig. la) with a staggered moment lying 
almost parallel to the Mn 3+ - O ap bond (where O ap is 
the apical oxygen in each pyramid). Neutron measure- 
ments show that the staggered moment is tilted away 
from the a-axis by about ±8 at low temperatures (14j . 
an effect which is presumably driven by a large single- 
ion anisotropy in the pyramid and strong AF exchange 
interactions along the chain (i.e., large J4 and Js>0 in 
Fig. la) [l|| . The effect of stacking AF zigzag chains 
along the 6-axis leads to a five member frustrated Mn 
spin loop in the afr-plane: Mn 4+ - Mn 3+ - Mn 3+ - Mn 4+ - 



Mn 3+ . Although the nearest-neighbor magnetic cou- 
pling in the loop is antiferromagnetic, the odd number of 
bonds leads frustrates the spins, preventing them from 
being be anti-parallel on every bond in the loop. Due to 
the relatively small exchange interaction between Mn 3+ - 
Mn 4+ ions (J3), the system forms a unique spin order- 
ing pattern [1J|: half of the Mn 3+ - Mn 4+ spin pairs 
across neighboring zigzags are approximately antiparal- 
lel, whereas the other half are almost parallel. Under this 
spin arrangement, exchange-striction between the Mn 3+ - 
Mn 4+ spin pairs shifts ions (mostly Mn 3+ inside pyra- 
mids) in a way that optimizes the spin-exchange energy: 
ions with antiparallel spins are pulled towards each other 
(green dashed circle in Fig. 1 a) , while those with parallel 
spins move apart (red dashed circle) [l8[ . This exchange- 
striction mechanism then results in canted staggered elec- 
tric dipoles that are aligned parallel to the Mn 3+ - O ap 
bond in each bi-pyramid, leading to an electric polariza- 
tion P along the 6-axis (Fig. 16). 

The application of high magnetic fields along the a- 
axis reveals unexpected magnetic and electric anomalies 
(Fig. 2); a sharp symmetric peak in the dM/dH curve 
(Fig. 2e) is accompanied by a sharp increase in e(H ) and 
a polarization P(H) which passes through zero (Fig. 2d 
and 2/). The critical fields for these anomalies increase 
with temperature, and more or less coincide. Similarly, 
the temperature of the e(T) peaks correspond to the tem- 
perature where the polarization P(T)=0 (Figs. 26 and 
2c). A contour intensity plot of the P(H) curves (Fig. 
le) confirms that the trajectory of the P(H)=0 (blue 
solid line) closely matches that of the P(T)=0 (red solid 
line), and the peak positions of e(H), dM/dH, and e(T) 
curves below ~20 K. 

We attribute the origin of the magnetic anomaly to a 
spin-flop (SF) transition in the AF zigzag spin chains. 
This is supported by an observation of a small M value 
and a zero-crossing of the linearly extrapolated line at 
H>H C , indicating the survival of the AF interaction be- 
yond H c (Fig. 2e, inset). While the SF is often man- 
ifested as a lst-order transition, in this case the M (H) 
curve does not show a discontinuity down to 0.5 K, as 
evidenced by a lack of divergence in the dM/dH curve. 
This characteristic is also observed in the electric prop- 
erties; the P(H) curve at 0.66 K does not show a dis- 
continuity and the e(H) curve does not show divergence 
even at T=0.66 K except sharpening on cooling. Fur- 
thermore, the e(H) curves measured in very slow H- or 
T-sweeps did not show any obvious hysteresis (See, a 
supplementary figure, Fig. SI). The lack of a sharp 
lst-order transition feature becomes more evident in the 
P(T) curves (Fig. 2c), in which P gradually changes 
its magnitude and sign on cooling. In addition, our pre- 
liminary T-dependent specific heat measurements across 
the P=0 line show no peak structure or discontinuity. 
These observations strongly suggest that the SF tran- 
sition and concurrent electrical transition occur via a 
smooth crossover down to temperatures as low as 0.5 K, 
during which the spin configuration changes continuously 
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FIG. 2: Electric and magnetic properties of BiM^Ov (a) 
H measured along the three crystallographic directions, (b) 
e(T) data down to 4 K at various magnetic fields. The 
e(H=0) peak near 40 K represents the FE transition, (c) 
P(T) curves at various magnetic fields, showing a crossover 
from positive to negative, (d) e(H) at fixed temperatures, (e) 
Field-dependent dM/dH curves at fixed temperatures. The 
inset shows the M(H) curve at T=0.5 K. (/) P(H) at fixed 
temperatures. The dashed line shows the calculated P curve 
as predicted by the spin-rotation model described in the text 
and in Fig. l(b-d). All dotted lines are guides for the eye. 



from the AF- (Fig. 16) to the SF-type (Fig. Id). 

In fact, continuous spin-flop transitions are well known 
in antifcrromagnets when the field H is applied at an an- 
gle to the easy axis: typically, a spin-flip crossover occurs 
when the angle between the field and the staggered mo- 
ment exceeds a certain critical value [l^, [2(|. This case 
has been well studied in uniaxial antiferromagnets such 
as GdA10 3 [H| and MnF 2 In the former, the critical 
angle starts from zero at Tn=3.2 K and grows as the tem- 
perature is lowered, reaching 10° at zero temperature. So 
long as the magnetic field inclination exceeds the critical 
angle, a cross-over is observed, but when the critical an- 
gle reaches the field-inclination angle, the system passes 
through a critical end point, located at a temperature 
T* in the H-T plane and at lower temperatures the SF 
sharpens into a first order transition. In spin flop tran- 
sitions, the location of T* is determined as a function 
of the initial angle between the spin easy axis and the 
applied field. 

We can easily generalize this understanding to the case 
ofBiMn205. Fig. la shows that in BiM^Os, the zig-zag 



chain of Mn spins (red arrows) subtend an angle of ±8 ° 
to the a-axis, providing an intrinsic "tilt" with respect to 
a field applied along the x axis. In this way, the zig-zag 
chain provides a natural mechanism for a spin-flop cross- 
over down to some critical temperature T* . This leads 
us to identify the P=0 trajectory in the phase diagram 
of BiM^Os (Fig. le) with the spin-flop crossover line 
above T>T*. Our data enable us to estimate the critical 
temperature T* is at, or below 0.5K. 

One of key findings of this work is the SF crossover 
can drive the continuous //-induced P reversal at T>T* 
via the exchange-striction mechanism. To understand 
the mechanism, consider the effect of a magnetic field 
H on a single AF zigzag chain. If we assume that the 
interchain coupling ( J3) can be neglected, then the dom- 
inant intrachain coupling (J4 and J5) will give rise to an 
almost antiparallel spin alignment. The magnetic free 
energy of the AF zigzag chain at zero temperature will 
then be governed by that of a pair of Mn 3+ spins in the 
bi-pyramid, consisting of the energies of the AF spin sub- 
lattice and the single-ion anisotropy of the Mn 3+ spins 
U : F M =-/i {x|| + {X± ~ X\\)sin 2 e}H 2 -Acos 2 (e - O ). 
Here, \\\ an d X-L are magnetic susceptibility values par- 
allel and perpendicular to the staggered moment, respec- 
tively. The parameters 9, 9o, and A denote the angle be- 
tween the Mn 3+ spin direction and the a-axis, the value 
of 9 at zero field and the single-ion anisotropy energy, 
respectively. The 9 value at each H leading to the min- 
imum energy can be obtained from 8Fm/ 'd9=0 and the 
resultant following relationship: 
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The above equation predicts that when 9 >0 (<0), 8 in- 
creases (decreases) continuously with H to cross the an- 
gle 45 (-45 °) at H c =y / Acos26 Q /ii (x_L ~ X||)- In otner 
words, the spin rotation occurs in a way that each AF 
zigzag chain rotates counter-clockwise or clockwise alter- 
natively to cross 45 or -45 at H c (dark blue arrows 
in Fig. la). When we consider the special ±45 con- 
figuration (Fig. lc), in which the angle between Mn 3+ 
and Mn 4+ spins becomes 90°, the exchange-striction is 
balanced and gives P—0. On further increasing H, \6\ 
increases to make the spin alignment almost perpendic- 
ular to the field (blue arrows in Fig. la). As a result, 
the relative orientation of neighboring Mn spins reverses 
as a function of field (see Fig. 1. b and d) giving rise to 
the sign reversal in the electric polarization P. 

For comparison with the experimental data, we have 
calculated predicted value of PocS 4+ -S' 3+ ~cos2# with 9 
determined in Eq. {l}. Here, S i+ and are the spin 
vectors of neighboring Mn ions between the AF zigzag 
chains. The dashed line in Fig. 2/ shows the calcu- 
lated P(H) curve with ^ #c=18.04 T and 9 =8 °, show- 
ing a qualitative agreement with the curve at T=0.66 K. 
Therefore, the rather continuous magnetic field-induced 
P reversal observed in BiMn 2 Os can be accounted for 
by the variation of exchange-striction degree driven by 
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the Mn spin-flop crossover. We further note that only 
BiMn 2 05 has a homogeneous positive P below Tq un- 
der small electric poling; all the other FM^Os develop 
additional FE domains with negative polarity and an- 
tiphase domain wall boundaries below Tic |22|. In turn, 
the above P reversal model can be applied to BiM^Os 
only, and it appears to explain why the experimental find- 
ing of P reversal as well as e increase near H c are unique 
to BiMn 2 5 [H. 

The magnetically induced polarization reversal we ob- 
serve is closely analogous to the field-induced metam- 
agnetic transitions observed in spin systems P, [24| and 
is most naturally regarded as a magnetic-field induced 
"metaelectric" transition [f| ■ As the inversion symmetry 
is already broken in both sides of the H c (T) boundary in 
Fig. 1/, the metaelectric transition involves only a change 
in the magnitude of the ferroelectric order parameter P. 
In this case, due to the coupling of P to the Mn spin mo- 
ment, the metaelectric transition occurs as a crossover 
down to very low temperatures (T=0.5 K). The system- 
atically increased sharpening observed in e{H) 1 dM/dH 
and P{H) curves upon cooling all indicate that T* is 
close to 0.5 K. If so, we would expect to observe crit- 
ical fluctuations associated with the critical end point, 
with a power-law evolution of the physical quantities in 
close proximity to T*. Our experiment confirms this ex- 
pectation as shown in Figs. 3a and 3b, where (e(H)- 
e(-ff=0))~ 1 vs. /i 2 / 3 and \P\ vs. h 1 / 3, plots produce quite 
linear behavior. Here h = \{H-H C / H c \ is the normalized 
deviation from the critical field F c =18.04 T. The linear- 
ity observed over a broad range of h (\h\=0. 004-0. 45 for 
e) means that a single power-law of Poc/i 1//3 can robustly 
explain the experimental data, even for H very close to 
H c . Furthermore, our attempts to fit the data at tem- 
peratures above 1.5 K reduce the linearity range for both 
e and P fits so that the power-law becomes most evident 
at the lowest temperature available, T=0.66 K, in our 
experiment. 

We find that the peculiar power-law of P and e near T* 
can be effectively described by a Landau-Ginzburg-type 
free energy F—Fme+Fe, where 

F ME = a(S 4 + • S 3 +)P = a(H - H C )P (2) 

and 

h ii 
F E = -(T-T*) n P 2 + -P 4 . (3) 
2 4 

The magnetoelectric free energy Fme consists of a cou- 
pling term between the spin orientations and the electric 
polarization, while the ferroelectric free energy Fe effec- 
tively describes the evolution of P as a crossover at T>T* 
and as a lst-order transition at T<T* under the assump- 
tion that the coefficient of P 2 changes its sign across 
T* . Here, n refers to an arbitrary exponent to describe 
temperature-dependence of physical quantities near T* 
and n—1 in the classical case. The second-derivative 
of Fe with respect to P determines the dielectric sus- 
ceptibility XcHP^HT-^Y'+ZuP 2 , where £=Xo (P) + l. 
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FIG. 3: Field-induced power-law evolution and e anomaly 
at very low temperature, (a) (e(H)-e(H=0))~ 1 vs. h 2 ^ 3 , 
where h=(H-H c ) / H c and fioH c = 18.04 T. A very slow ramp- 
ing rate (0.1 T/min) was used for the e(H) measurement. 
The red dashed lines are linear guides passing through the 
origin. The linear relationship was valid in a wide |ft|=0.004- 
0.45 (i.e., ^ |#-#c|=0.1-8 T). (b) Plot of \P\ vs. h 1/A . The 
linear regime of h was \h\=Q. 02-0.45 (^ o jF-F c |=0.3-8 T). To 
increase the signal-to-noise ratio, it was necessary to choose 
a rather fast ramping rate (5 T/min) for the P(H) measure- 
ment, which made the P data less reliable in a small h region, 
(c) Detailed e(T) curves down to T~0.6 K at fixed H near 
H c . The shape of the e(T) curves above the maximum tem- 
perature was similar to that predicted by Barrett's theory 
(Fig. S2). 



The sharp increase in Xc observed in experiments can 
be understood as the closeness of the system to the 
critical end point T* . In our Landau-Ginzburg model, 
when the temperature is approximately T*, the physics 
is driven by the interplay between the dominant quar- 
tic P 4 term and linear P term. Under this assumption, 
d(FuE+F E )/dP\ E=0 =0 provides Poc( J ff- J ff c ) 1 / 3 . Fur- 
thermore, d 2 (F M E+F E )/dP 2 =x~ 1 =3uP 2 (x(H-H c ) 2 / 3 . 
These predicted exponents are consistent with those de- 
termined experimentally (Figs. 3a and 36), clearly sup- 
porting the existence of T* near 0.66 K. 

The low value of T* found in our experiments suggests 
that BiMn205 lies close to a a metaelectric quantum 
critical point, where T* is actually zero. Closely anal- 
ogous behavior is well-known in the case of the itinerant 
metamagnet Sr3Ru20y [l|. In this respect, several novel 
features of the present data particularly prompt further 
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notice. First, the observation of the robust power-law 
of Poc/1 1 / 3 suggests a unique mean-field type while the 
system is expected to have a 3D Ising universality class 
(Pocft. 1 / 4 ' 82 ) near the classical critical end point as seen in 
the liquid-gas transition. Moreover, the trajectory of the 
crossover temperature near the zero temperature limit 
following (H-H c ) 1 / 2 (Fig. le) again constitutes another 
mean-field type behavior. Thus, all these observations 
suggest that the effective dimension of the system has 
been lifted up to 4 or higher so as not to produce the 
expected 3D fluctuation. One possible origin for such a 
dimensional increase is in fact the dynamic exponent ef- 
fect expected due to the prevailing quantum critical end 
point fluctuation. 

The insights gained by our measurements and analy- 
sis above suggest that the multiferroic BiM^Os provides 
a rare opportunity to find a multi-ferroic quantum crit- 
ical point and to test its influence on ferroelectric ma- 
terials. Supporting this novel conclusion, detailed e(T) 
shape realized near H c (Fig. 3c) further reveals a clue 
toward detecting the quantum critical signature. When 
H approaches 18.04 T from below, systematic increase 
of e(T=0.66 K) occurs. In particular, at /x i? c =18.04 
T, e continues to increase as T decreases from 20 K and 
shows a broad maximum around 4 K, followed by a slight 
decrease of 0.5 % down to ^0.6 K. The behavior qual- 
itatively resembles that of the s(T) curve in quantum 
paraelectric systems, such as SrTiOs [If|[2(|[27, 28]. Fur- 
thermore, the increase in e(T) between 4 and 20 K is not 
proportional to the inverse of temperature, as expected 
in classical paraelectrics. Instead, it is closer to the shape 
predicted by Barrett's theory (Fig. S2) [29| . indicat- 
ing that the enhanced lattice fluctuation in close vicin- 
ity to T* is linked to the quantum mechanical population 
of the relevant phonon mode. Therefore, the physics of 
BiMn2 05 appears to involve quantum zero point lattice 
motion at low temperatures as does SrTiC>3. The unique- 
ness of BiM^Os lies in the magnetic-field-tunability of 
its quantum-paraelectric properties. 

Our observations raises the interesting possibility that 
instabilities into new quantum states of matter might 
develop in close vicinity to the multi-ferroic critical end 
point, as they do close to the magnetic quantum critical 
end point in S^R^Oy. In future work we plan to in- 
vestigate this possibility by fine-tuning the T* position, 
either by by rotating the field direction or by applying 
external pressure. 



Methods 

Single ^crystal samples were grown using the flux 
method Q. All the electrical properties were measured 
along the 6-axis while the magnetic field was applied 
along the a-axis up to 35 T at the National High Mag- 
netic Field Laboratory (NHMFL). Dielectric constants 
were measured by a capacitance bridge (GR1615A or 
AH2500A) at 1 kHz [23[. P was obtained by integration 
of the pyro- and magnetoelectric currents measured by 
use of an electrometer (Keithley 617). Given the limited 
magnet time for high-field measurements, the continuous 
application of a small electric field bias of E=2 kV/cm 
was quite useful in quickly determining the P evolution. 
In particular, the value of E=2 kV/cm was sufficient to 
fully polarize the sample on cooling through Tq but small 
enough to produce a P close to its spontaneous value. In 
several field and temperature sweeps, it was verified that 
the spontaneous P measured under zero electric field af- 
ter poling with i?=2 kV/cm is close to the P measured 
under the same bias. A precalibrated capacitance sensor 
was used for accurate T-sweeps between 0.6 and 20 K 
near H c (Fig. 3c), while a Cernox sensor was mostly used 
for the T- and .ff -sweeps above 4 K. For H-swept mea- 
surements at T"=0.66 K and 1.5 K, temperatures were 
stabilized by controlling the vapor pressure. Specific heat 
and magnetocaloric effects were measured with a plastic 
calorimeter. Temperature-dependent DC magnetization 
was measured using a vibrating sample magnetometer 
in a Physical Property Measurement System (Quantum 
Design). A pickup-coil magnetometer was used to mea- 
sure the if-dependent magnetization in a pulse magnet 
at NHMFL-Los Alamos National Laboratory. 
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